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Single particle interference lies at the heart of quantum mechanics. The archetypal double-slit 
experiment^ has been repeated with electrons in vacuum^'^ up to the more massive Cgo molecules^. 
^-^ Mesoscopic rings threaded by a magnetic flux provide the solid-state analogous^'®. Intra- molecular 
interference has been recently discussed in molecular junctions^'^!i^ i^°i^^ . Here we propose to exploit 
interference to achieve all-electrical control of a single electron spin in quantum dots, a highly de- 
psj sirable property for spintronic a^^i^^i^^ and spin-qubit application o^^i^^'-'^'^'-'^^'-'^^ . The device consists of 
5_i an interference single electron transistor (ISET)-*^"'^^, where destructive interference between orbitally 
degenerate electronic states produces current blocking at specific bias voltages. We show that in the 
presence of parallel polarized ferromagnetic leads the interplay between interference and the exchange 
coupling on the system generates an effective energy renormalization yielding diff"erent blocking biases 
for majority and minority spins. Hence, by tuning the bias voltage full control over the spin of the 
,— I trapped electron is achieved. 

The all-electrical solutions to the challenge of single spin control that have been propose d^^i^^'^^i^^ and realize d^^i^° 
are based either on spin orbit-couplin g °' °i 'i °i^ or on tunneling-induced spin splitting in the Kondo regimes^. Our 
proposal relies on the current blocking occurring in an ISET due to interference between degenerate states. The 
C/3 , conditions for interference blocking are very generic^^ and admit several different realizations. We consider here for 
^ • clarity a benzene and a triple-dot ISET, Fig.[TJ Both are described by the Hamiltonian: 

^ H = Hgys + -ff loads + Ht, (1) 

where Hgys represents the central system and also contains the energy shift operated by a capacitively coupled gate 
r ■ electrode at the potential Vg. The Hamiltonian Hgys is in both cases invariant with respect to a discrete set of rotations 
■ O ■ around the vertical axis passing through the center of the system. This fact allows a classification of its eigenstates 
in terms of the z component of the angular momentum £ and also ensures the existence of degenerate states with 
different £. Then, a generic eigenstate is represented by the ket \NiaE) where N is the number of electrons on the 
system, a is the spin and E the energy of the state. When degenerate states participate to transport they interfere 
since, like the two paths of the double-slit experiment, they are occupied simultaneously by the travelling electron, 
but in different superpositions under diverse transport conditions, -fficad describes the ferromagnetic leads with equal 
(for simplicity) parallel polarization P and with a difference eVb between their electrochemical potentials. Finally, 
\^ I accounts for the weak tunnelling coupling between the system and the leads, characteristic of SETs, and we consider 
the timnelling events restricted to the atoms or to the dots closest to the corresponding lead (Fig. [T]). We explicitly 
' consider the Coulomb interaction only in the central part of the device (see the supplementary material Q]) due to 
, the strong confinement experienced there by the electrons while, apart from the polarization assumption, we assume 

■ a non interacting approximation for the leads. 

^\ I In the weak coupling regime the current essentially consists of sequential tunnelling events at the source and drain 

■ lead that increase or decrease by one the number of electrons on the system. The different panels of Figs. [5] and [3] 
^ ' show the current through the benzene and triple dot ISET, respectively, as a function of bias and gate voltage. As 

• I— I 1 in all SETs at low bias so called Coulomb diamonds, where transport is energetically forbidden, occur. Within the 
' diamonds the particle number is fixed as indicated in the figures. Only exceptions are the charge degenerate points 
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FIG. 1: Two examples of interference single electron transistors (ISETs): a benzene molecular junction contacted in the meta 
configuration (A) and a triple quantum dot artificial molecule (B). The source and drain are parallel polarized ferromagnetic 
leads. 
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where two diamonds meet. Here the energy difference of two ground states with consecutive particle numbers is 
equal to the equilibrium chemical potential of the leads. At finite bias the incoming electrons have enough energy to 
overcome the level spacing and the Coulomb repulsion and the current flows. As a signature of the new states that 
enter the bias window, by increasing the voltage the current typically increases steplike. 

In ISETs an exception to this picture is represented by the interference blockade where the current decreases for 
increasing bias generating negative differential conductance (NDC) and eventually vanishes (see green lines in the 
panels B and C of Fig. [5] and [3]) . Panels B in the same figures indicate moreover that, for a given gate voltage and in 
absence of polarization in the leads, the current is blocked only at one specific bias voltage. For parallel polarized leads, 
however, at a given gate voltage, the current is blocked at two specific bias voltages, one for each spin configuration 
(panels C). As demonstrated below, the blocking of the minority electrons occurs for the smaller bias voltages. As 
such full control of the spin configuration in the ISET can be electrically achieved. The interference blockade and its 
spin selectivity is also demonstrated in panels A and B of Fig. [H Along the dotted (dashed) line a majority (minority) 
spin electron is trapped into the molecule. The molecular spin state can thus be manipulated simply by adjusting the 
bias across the ISET. In the following we discuss the physics of the spin-selective interference blocking and present 
the necessary ingredients for its occurrence. 

This novel blocking is explained by the presence of an A^-particle non-degenerate state and two degenerate + 1- 
particle states that simultaneously contribute to transport. It also requires that the ratio between the transition 
amplitudes jai (* = 1, 2, a — L, R) between those N- and N + 1-particle states is different for tunneling at the left 
(L) and at the right (R) leadii: 



This condition is fulfilled in both cases presented in Fig. [T] due to the geometrical configuration of the left and right 
lead. Due to condition ([2]), the degenerate states interfere among themselves such to form pairs of blocking and 
non-blocking states. The blocking state is only coupled to the source lead (panels C and D of Fig. |3|) while the 
non-blocking one to both source and drain. An electron that populates the blocking state can neither leave towards 
the drain nor, at high enough bias, return to the source since all energetically available states are there filled. The 
non-blocking state is thus excluded from the dynamics and the current vanishes. 

As such we would conclude that the interference blocking is a threshold effect and the current remains blocked until 
a new excited state participates to the transport. However, as shown in Fig. [2] and [3l the current is blocked only at 
specific values of the bias voltage. The explanation of this phenomenon relies on two observations: i) The blocking 
state (Fig. |4]) must be antisymmetric with respect to the plane perpendicular to the system and passing through its 
center and the atom (quantum dot) closest to the drain; this state is thus also an eigenstate of the projection of 
the angular momentum in the direction of the drain lead^i. At positive (negative) bias voltages we call this state 
the i?(_L)-antisymmetric state \ipR(L),a)- ii) The coupling between the system and the leads not only generates the 
tunneling dynamics described so far, but also contributes to an internal dynamics of the system that leaves unchanged 
its particle number. In fact the equation of motion for the reduced density matrix p of the system can be cast, to 
lowest non vanishing order in the coupling to the leads, in the form: 



The commutator with -ffgys in El- (El represents the coherent evolution of the system in absence of the leads. The 
operator £tun, describes instead the sequential tunnelling processes and is defined in terms of the transition amplitudes 
7ai between the N and A'^ -I- 1 particle states like the ones introduced in equation ([2|) . Eventually Hcs renormalizes 
the coherent dynamics associated to the system Hamiltonian. It reads: 



where La is the projection of the angular momentum in the direction of the lead a and, for paramagnetic systems, 
it does not depend on the spin degree of freedom a. Moreover, tOaa is the frequency renormalization given to the 
states of spin a by their coupling to the a lead. Equation ^ is an example of Bloch-Redfield equation describing 
the dynamics of a system coupled to thermal baths. A more detailed version of ^ is presented in the supplementary 
material [2] 

For sake of simplicity we give in the following the explicit form of the transition amplitudes jai, of the operator 
La and of the associated frequency Waa- only for the benzene ISET and for the ground state transition 6g — > 7g that 
is characterized by interference blocking. The argumentation is nevertheless very general and can be repeated for all 
the systems exhibiting rotational symmetry. The transition amplitudes read: 
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FIG. 2: Benzene ISET: polarized vs. unpolarized configuration. Panel A - Current vs. bias and gate voltage for unpolarized 
leads. Panel D - Current vs. bias and gate voltage for polarized leads (polarization P — 0.85). Panels B and C - Blow up of the 
6-^7 particle transition for both configurations. The unpolarized case shows a single current blocking line and the trapped 
electron has either up or down polarization. The polarized case shows two current blocking lines, corresponding to the different 
spin of the trapped electron. The current is given in units of e/F where F is the bare average rate (supplementary material [4} , 
and the temperature ksT = O.Olfe where b is the hopping parameter (supplementary material nj) 



lai = (6g00|dM<T|7g^cr)e" 



(5) 



where \7g ia) are the orbitally degenerate 7 particle ground states, £ = ±2 the z projection of the angular momentum 
in units of H and dj^^ destroys an electron of spin tj in a reference carbon atom M placed in the middle between the 
two contact atoms. Moreover, 0q is the angle of which we have to rotate the molecule to bring the reference atom M 
into the position of the contact atom a. The present choice of the reference atom implies that (f>L — —4>r — f ■ 
the Hilbert space generated by the two- fold orbitally degenerate \lgia) the operator reads: 
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(6) 



For a derivation of ([6|) see the supplementary material [31 The frequency ujaa is defined in terms of transition 
amplitudes to all the states of neighbour particle numbers: 



•'au — 7^ i aa' 



(JMdiiME]mE]KiJlama)p^{E^~E) 



(7) 



where the compact notation \N{E}) indicates all possible states with particle number N and energy E, Pa{x) = 



^{x — fia) where /3 = l/k^T, T is the temperature and ip is the digamma function. Moreover T\ 
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FIG. 3: Triple dot ISET: polarized vs. unpolarized configuration. Panel A - Current vs. bias and gate voltage for unpolarized 
leads. Panel D - Current vs. bias and gate voltage for polarized leads (polarization P = 0.7). Panels B and C - Blow up of 
the 6^5 particle transition for both configurations. The selective spin blocking is analogous to the one of the benzene ISET 
(Fig.[2j. 

^\t\'^Daa' is the bare tunneling rate to the lead a of an electron of spin a', where t is the tunnelling amplitude 
and Daa' is density of states for electrons of spin a' in the lead a at the corresponding chemical potential /Xq. Due 
to the particular choice of the arbitrary phase of the 7 particle ground states, ujaa does not depend on the orbital 
quantum numbers £ and m. It depends instead on the bias and gate voltage through the energy of the 6, 7-ground 
and 8 particle states. In Fig. O the black curve depicts lulo- as a function of the bias in absence of polarization: the 
frequencies corresponding to the two spin species coincide and thus vanish at the same bias. The same condition, 



UJLa = 0, (8) 

also determines the bias at which the current is completely blocked. In fact, at that bias the effective Hamiltonian 
contains only the projection of the angular momentum in the direction of the right lead (the drain) and the density 
matrix corresponding to the full occupation of the 7 particle i?-antisymmetric state (p = {ipR^ a)(V'-R, a|) is the stationary 
solution of Eq. As we leave the blocking bias the effective Hamiltonian contains also the projection of the angular 
momentum in the direction of the left lead and the i?-antisymmetric state is no longer an eigenstate of H^s- The 
corresponding density matrix is not a stationary solution of ^ and current flows through the system. The L ^ R 
symmetry of the system implies, for negative biases, the blocking condition Wflo- = 0. 

All-electric-spin control is achieved, in an ISET, only in presence of ferromagnetic leads and with exchange interac- 
tion on the system. By manipulating ([7]) it is possible to show that the frequency splitting LUa^ — tOai is proportional 
to the polarization in the a lead, but vanishes in the absence of exchange interaction on the system capable to lift 
the singlet-triplet degeneracy of the excited 6 and 8 particle states (see the supplementary material [4)1 . In Fig. [5] we 
show the frequencies ojl^ — vs. bias voltage also for a finite values of the polarization P calculated for the benzene 
ISET, where exchange splitting is ensured by the strong Coulomb interaction on the system. The interference blocking 
conditions t^Lo- — for the L ^ R current are satisfied at different biases for the different spin species. The dotted 
and dashed lines in Fig. |4] are the representation of the relations ujli = 0,^Li = as a function of the bias and 
polarization, respectively. 

In previous studie&i^ we have shown that the interference current blocking does not depend on the perfect symmetry 
of the system but rather relies on the existence of quasi-degenerate states in which the energy splitting is smaller than 



5 



Current Minority blocl<ing 




0.4 0.8 1.2 1.6 

Vb [|b|/e] 



FIG. 4: Spin control. Panel A - Current through the benzene ISET vs bias and polarization at the 6 — > 7 electrons transition. 
Panel B - Population of the majority spin 7 particle state. The two zero current lines at high bias correspond to the maximum 
or minimum population of the 7 particle majority spin state and thus identify the spin state of the trapped electron on the 
molecule. Panels C and D - Schematic representation of the spin selective blocking corresponding to the dashed (C) and dotted 
(D) lines of the panels A and B. 
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FIG. 5: Blocking condition. Renormalization frequencies lulo- of a benzene ISET as function of the bias and for different lead 
polarizations. The current blocking condition ujLcr = is fulfilled at different biases for the different spin states. 



the tunnelling coupling to the source and drain leads. In the proposed structures the degeneracy is associated with 
the rotational symmetry and it has the advantage of a simple geometrical realization of the interference conditions 
^ . Nevertheless the effect is more general and any other structure exhibiting orbital degeneracy is a good candidate 
for an ISET. 
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SUPPLEMENTARY MATERIAL 
1. The system Hamiltonian 

The Hamiltonian that describes both systems represented in Fig. [Dreads 

i(T ia 

+UJ2{n^1 - {n^l- ^) (9) 

i 

+y ^ + n^i - 1) (n,+ii- + n^+n - 1) , 

i 

where dt^ creates an electron of spin cr in the orbital of carbon i or in the ground state of the quantum dot i and 
i = 1, . . . , 6(3) runs over the six carbon atoms (three quantum dots) of the system. Moreover riia = d\^d^^. The 
effect of the gate is included as a renormalization of the on-site energy ^ = — eVg {Vg is the gate voltage) and 
we conventionally set Vg = at the charge neutrality point. The parameters that we have used are b = 2.5eV, U = 
9eV, V = 6eV. 



2. The generaHzed master equation 

We describe the dynamics of the system with a generalized master equation (GME) for the reduced density matrix 
p. This equation is obtained from the Liouville equation for the full density matrix as a perturbation to the lowest 
non vanishing order in the coupling to the leads by tracing out the leads degrees of freedom. A generic formulation of 
the GME for a SET in presence of degeneracies or quasi-degeneracies can be found elsewhere e.g}^. We concentrate 
here on the range of gate and bias voltages at which the dynamics is restricted to transitions involving the |6g00) and 
\7g£a) many particle states of the benzene ISET. 

The seven particle states are spin and orbital degenerate. The general theory of the GME would require a priori 
to keep thus a full 4x4 density matrix describing the 7 particle subspace. In presence of parallel polarized leads, 
though, the coherences between different spin degrees of freedom can be neglected since spin is always conserved 
by the electrons while travelling through the device. The GME can thus be written in terms of the nine variables 
collected in the 1x1 matrix p^o and the two 2x2 matrices p^""^ with a |. Due to the rotational symmetry of the 
system it is more convenient to refer to another set of variables, namely to describe the dynamics in terms of the 
occupation probabilities Wg, Wro- and the expectation values of the different projections of the angular momentum 
for the system. The new set of variables is: 

W6 = p^^ 

W7a = Trjp^''^} 
= Tr{L,p7«"} 
= Tr{Ly^^}. 

The operator Lz is the generator of the set of discrete rotations around the axis perpendicular to the plane of the 
benzene molecule that bring the molecule into itself and can be written within the 7 particle Hilbert space spanned 
by the vectors \7g£(7) as = —h\£\az, where az is the third Pauli matrix. The operator La generates, in the same 
space, the discrete rotations around the axis in the molecular plane and passing through the center and the atom 
closest to the contact a. Finally, the dynamics for the variables introduced in Eq. (|10p is given by the equations: 
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We =2^r„, [-f+iAE)We + /-(Ai?)L„,] , 

acr 

W7. =2^r„, [f+iAE)We - /-(A£;)L„,] , 

a 

= - 2r,,/-(A£;) + 2 {r„,/+(Ai;) + ra,/+(As) cos2[|^|((/,„ - ^^)]} We 



(11) 



+ T^^f^iAE) sin'micj,^ - (/.a)]W7. - r^./^ (Aii;)(L„, + L^.) + ^^^^M^^^^c^a.W, 
W = - ^ r„,/-(A£;)L,, - 2tan[|^|((/.L - 0fl)](wLa - ujr„){W7, - Ll^ - Lr,) 

a 

- 2c0t[\e\{(f>L ~ <t)R)]{^^La + UJB,a)(LLa ~ ^Ra), 

where V^a — TqctI (^gOOMacrl^g^cr) p is the tunneUing rate at the lead a involving the ground states with 6 and 7 
particles. Terms describing sequential tunnelling from and to the lead a are proportional to the Fermi functions 
/+(a;) := ](x - fia) and f~{x) := 1 - respectively, and AE = E^g - E-jg + eVg where Ef^g and Ejg are the 

energies of the 6 and 7 particle ground states. Finally with a we mean the lead opposite to the lead a. By using the 
expression \l\ (to be substituted with 2 for the 6 — > 7 particle transition) we maintained the generality of the equations. 
The replacement |-^| = 2 — > 1 and the appropriate redefinition of AE is enough to treat the 6^5 transition. Another 
important generalization concerns the position of the leads. The para (0^ — (j)R = tt) and ortho {(pL — (t>R ~ '^1'^) 
configuration are also treated within the same equations. In particular one can see that all the terms containing the 
renormalization frequencies drop from the equations in the para configuration and that the equations for the ortho 
and meta configuration coincide. 

3. Matrix form of the operator La 

The explicit form of is given in Eq. We give here its derivation. It is convenient, for this purpose, to 
choose the arbitrary phases of the states \7gia) in such a way that the rotation of tt around the axis passing through 
a reference atom M and the center of the molecule transforms \7g£a) into — |7g — £(7). In other terms 

exp(z7r^^) ^ (12) 

where ax is the first Pauli matrix. The relation is in fact an equation for Lm and the solution reads: 

Lm^^(]]). (13) 



2 V 1 1 , 

Eventually we obtain La by rotation of Lm in the molecular plane, namely: 



La = e-i^=^-LMei^°^^ = | ^ , ) , (14) 



where is the angle of which we have to rotate the molecule to bring the reference atom M into the position of the 
contact atom a. 
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4. The spin splitting of the renormalization frequencies 

The spin splitting of the renormaUzation frequencies is obtained from Eq. ([7]) . By introducing the average bare rate 
r = °^ ^ , for simplicity equal in both leads, and using the fact that benzene is paramagnetic we get: 

TT ^ — ' 

{E} 

'{7gl T \d,,^\8{E}){8{E}\dl^\7gm ^)p^{E ~ ErJ 

+ {7gl T \dl^\6{E}){6{E}\d,,^]7gm T)Pa(i?7,- E) ^^^^ 
-{7g£ T \d,,^\8{E}){8{E}\d\,^\7gm ])p^{E - EjJ 
-(7/ T \d\,^\6{E}){6{E}\dj,,^\7gm ^)p^{Ej- E)] , 

where one appreciates the linear dependence of the spin splitting on the lead polarization P^. The first and the third 
term of the sum would cancel each other if the energy of the singlet and triplet 8 particle states would coincide. 
An analogous condition, but this time on the 6 particle states, concerns the second and the fourth terms. For this 
reason the exchange interaction on the system is a necessary condition to obtain spin splitting of the renormalization 
frequencies and thus the full all-electric spin control. 



